Renormalization Group Study of Magnetic Catalysis in the 3d Gross-Neveu Model 
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Magnetic catalysis describes the enhancement of symmetry breaking quantum fluctuations in 
chirally symmetric quantum field theories by the coupling of fermionic degrees of freedom to a 
magnetic background configuration. We use the functional renormalization group to investigate this 
phenomenon for interacting Dirac fermions propagating in (2+l)-dimensional space-time, described 
by the Gross-Neveu model. We identify pointlike operators up to quartic fermionic terms that can 
be generated in the renormalization group fiow by the presence of an external magnetic field. We 
employ the beta function for the fermionic coupling to quantitatively analyze the field dependence of 
the induced spectral gap. Within our pointlike truncation, the renormalization group flow provides 
a simple picture for magnetic catalysis. 



I. INTRODUCTION 

In the past 20 years it has become evident that in cer- 
tain planar condensed matter model systems defined on 
some lattice, the low-energy physics can efficiently be de- 
scribed in terms of Dirac fermions. These "relativistic" 
degrees of freedom typically arise by linearizing a given 
dispersion relation around a special set of points in the 
Brillouin zone, where either a linear band crossing is re- 
alized or the single-particle gap function vanishes at iso- 
lated points across a given Fermi surface. 

The most prominent and well-investigated example 
is probably the Hubbard model on the honeycomb lat- 
tice This model describing the hopping of fermions 
on a hexagonal lattice with possibly strong on-site in- 
teractions serves as a minimal model for the electronic 
correlations in monolayer graphene and similar systems. 
Once the naive continuum limit has been taken, one can 
expect a faithful representation of the low-energy sec- 
tor by a quantum field theory for Dirac fermions moving 
in (2 4- l)-dimensional space-time in the proximity to a 
second order quantum phase transition. The lattice sym- 
metries of the Hubbard model get replaced by an "emer- 
gent" Lorcntz symmetry and the free part of the Dirac 
action is invariant under continuous chiral transforma- 
tions. Note that for the chiral symmetry to be realized, 
a reducible four-component representation of the Dirac 
algebra is necessary, which can be formed from the direct 
sum of two inequivalent two-component representations. 
In the case of an irreducible representation, the chiral 
symmetry effectively corresponds to part of the flavor 
symmetry. In the following, wc will consider the chirally 
invariant case. This asymptotic chiral symmetry is also 
realized for example in d-wave superconductors 0, 0] or 
quantum spin-liquid phases 0,0. In both cases, how- 
ever, a symmetry of the underlying parent state is al- 
ready assumed to be broken, in contrast to the case of 
the Hubbard model on the honeycomb lattice. Another 
situation where no symmetry breaking is necessary are 
topological insulators and superconductors, which sup- 
port Dirac fermions as surface states 0, Q • 

Once interactions are included in such a low-energy ef- 
fective quantum field theory, we expect that the chiral 



symmetry is spontaneously broken as soon as the associ- 
ated coupling exceeds a critical value. The fluctuation- 
driven interactions compatible with the symmetries de- 
termine in which channel the formation of a conden- 
sate, i.e., a finite expectation value of an order-parameter 
field, can take place. In this paper, the focus of our 
investigation is on the Gross-Neveu model with a re- 
ducible representation of the Dirac algebra and a purely 
pointlike scalar-scalar interaction. This theory realizes 
a system with discrete chiral symmetry. Quantum phase 
transitions falling into the Gross-Neveu universality class 
for two fermion flavors Nf = 2 are, for example the 
excitonic Q and antiferromagnetic instabilities (lo| in 
graphite and graphene, respectively. The universal prop- 
erties of secondary d- to d + is-wave pairing transitions 
in nodal d-wave superconductors are also expected to be 
described by the Gross-Neveu model 0, 01 ■ Chiral sym- 
metry breaking is known to occur for all flavor numbers 



A^f whereas in models with continuous chiral sym- 



metry and interactions in a vector-vector channel such 
as the Thirring model a chiral condensate is stable 
only up to a critical flavor number. 

Remarkably, when interacting planar Dirac fermions 
are subject to an external perpendicular magnetic field, 
symmetry breaking accompanied by the generation of a 
mass gap for the fermionic degrees of freedom occurs for 
all values of the interaction strength [l^. Indeed, this 
effect known by the name of magnetic catalysis is not 
particular to models defined in three space-time dimen- 
sions, but has also been investigated in four dimensional 
quantum field theories However, what is special 

about three space-time dimensions is the fact that even 
in the absence of interactions a finite value of the chiral 
condensate is generated by applying an external magnetic 
field [l^ . So also for a free system the (continuous) chiral 
symmetry is spontaneously broken. This phenomenon is 
sometimes referred to as a 'quantum anomaly' [l6| . since 
it unexpectedly leads to the generation of a non- vanishing 
expectation value for a certain operator. For condensed 
matter model systems, possible applications of this ef- 
fect have been studied for a variety of low-dimensional 
physical realizations, ranging from magnetic field induced 
insulating behavior in semi-metallic structures [itI [isj 
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to anomalous Hall plateaus in graphene [16|, Il9l - l23ll and 
transport properties in c?-wave superconductors j24M2(: 
It is interesting to note that magnetic catalysis can also 
go along with an induced anomalous magnetic moment 
of the fermionic excitations 
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The interplay between strong magnetic fields and 
chiral-symmctry breaking has also become of topical rel- 
evance in particle physics. For instance, non-central 
heavy-ion collisions go along with extreme magnetic fields 
[28l | that may even interfere with the topological struc- 
ture of the QCD vacuum [2^, [s^] or simply take a 
parametric influence on the chiral phases of QCD (sTI - 
1371 ■ Even though magnetic fields in strongly- interacting 
fermionic systems generically support chiral-condensate 
formation at least in mean-field-type approximations, 
further interactions such as gluonic back-reactions may 
also lead to inverse effects j38| . 

There exist already a few works in the literature on 
the large- TVf results for the Gross- Neveu model j39l - l4]| . 
Beyond leading-order results in a l/A''f-expansion, mag- 
netic catalysis was studied both with truncated Dyson- 
Schwinger equations (in d = 3 [l3, 113 and d = 4 [TS. lisl. 

and Wilson-style renormalization group (RG) equa- 
tions (in c? = 4 (45l - [47| ) in the context of fermionic mod- 
els with short-ranged interactions and QED. However, 
a clear RG picture in terms of a fixed-point analysis of 
the beta function appears to be still missing. Moreover, 
the inclusion of purely magnetically induced operators in 
the RG flow and their effect on the running coupling has 
not been considered so far. In this paper, we provide 
a clear renormalization mechanism for magnetic cataly- 
sis in the Gross-Nevcu model. Similar observations are 
made in an RG study of magnetic catalysis in a QCD 
low-energy model [i^. Moreover, we identify terms that 
are compatible with the symmetries of the Gross-Neveu 
action in the presence of an external magnetic field. In 
principle, all operators compatible with the symmetries 
of the initial action can be generated under RG transfor- 
mations. Therefore, our analysis addresses the question 
as to whether the operator content of fermionic effective 
actions is modified in the IR by the very presence of such 
gauge backgrounds. 

The paper is organized as follows. In Sect.|TTl we define 
the three-dimensional, chirally symmetric Gross-Neveu 
model in an external magnetic field. In Sect. IIIIl we col- 
lect the transformation rules under both discrete symme- 
try and chiral transformations, which we will then em- 
ploy in Sect. HV] to study the Gross-Neveu theory space. 
Put differently, we give the transformation properties of 
all bilinear operators and list all quartic operators com- 
patible with the Gross-Neveu symmetries. This includes 
the set of magnetically induced operators (up to a specific 
mass dimension) , the presence of which is completely sus- 
tained by the external magnetic field. To keep the pre- 
sentation self-contained, in Sect. |V] we recapitulate the 
functional RG essentials for the one-particle irreducible 
(IPI) scheme. In Sect. I VII we discuss the effect of the 
magnetic field on the behavior of the beta function for 



the fermionic coupling in the pointlike approximation in a 
qualitative manner. A quantitative discussion of the field 
dependence of the dynamically generated fermion mass 
inferred from the beta function is deferred to Sect. I VIII 
We conclude with Sect. IVIIII where we summarize our 
findings and give an outlook to our future work. Tech- 
nical details and the derivation of the beta function are 
provided in Appendices EUCl 



II. GROSS-NEVEU MODEL IN AN EXTERNAL 
MAGNETIC FIELD 

We consider the Gross-Nevcu model [i^ within func- 
tional integral quantization. The microscopic degrees of 
freedom are represented by Grassmann-valued fields de- 
fined over three-dimensional Euclidean space-time. For 
studying the effects of magnetic catalysis, we mini- 
mally couple the Dirac fermions to a background gauge- 
potential that realizes the external magnetic field. The 
defining local action to be quantized is given by 

S[tpj,ipj,A] ^ / C{Tjjj,ipj,d^,'ipj,d^'ipj,A) (1) 



with the Lagrangian density 



(2) 



where — J df^x is a shorthand for the integral over the 
d-dimensional Euclidean space-time. Here, 



Ip[A] = 7m (^m - 19-4^(2^)) , M - 0, 1, 2, 



(3) 



denotes the covariant derivative acting on the spinor 
fields with q the respective charge under the Abelian U(l) 
gauge group. We consider A^f different flavor species in- 
dexed by j = 1, . . . , A'f , but keep the flavor number arbi- 
trary. This model enjoys a global \]{Ni) flavor symmetry. 
The canonical mass dimensions of the fields are given by 
[■0] = and [A\ = This renders the gauge charge 



a dimensionful quantity with [q] 



As mentioned 



in Sect. U we use a reducible four-component represen- 
tation for the gamma matrices in this work, i.e., = 4 
where = trl4 denotes the dimension of the represen- 
tation space of the Dirac algebra and tr is the trace over 
spinor indices. Our choice for the 4x4 representation of 
the Dirac algebra is the so-called chiral one. Here, the 
75 matrix is in block-diagonal form. The matrices are 
explicitly given by 



70 = ^^2 <X) T3 , 71 = T2 (X) Ti , 72 = T2 (Xl T2 



and 



73 = n (Xi To , 75 = T-3 (K) To . 



735 = 17375 



(4) 



(5) 



Here, the {Ti}'s denote the Pauli matrices which satisfy 
TiTj = SijTQ + ietjkTk, with i,j,k = 1,2,3 and To = I2 



3 



is a 2 X 2 unit matrix. The gamma matrices satisfy the 
anti-commutation relations 

{7p;7i^} = 2(5^^.14, (6) 

where ^,,v = 0,1,2. The matrices 73 and 75 anti- 
commute with all 7^, ^ = 0, 1, 2. In this representation, 
72, 73 and 75 are symmetric and real, while 70, 71 and 735 
are antisymmetric and purely imaginary 4x4 matrices. 



One easily verifies the following identities 

= (5^1/735 + '-^^^ivala, (7) 

'^Ibliilv = 1(5^,. 75 + iep,.^737fT, (8) 

i737A<7>^ = Wp,.73 - '^<^iivalbla, (9) 



where the second and third lines follow from the first 
one by left or right multiplication with an appropriate 
gamma matrix and making use of the Dirac algebra, and 
e^ii/o- is the completely antisymmetric tensor. 
We chose the vector potential in the gauge 

A{x) ^ {Q,Q,xiBf . (10) 

Within the three-dimensional formulation, the magnetic 
field appears as a pseudo-scalar quantity. This corre- 
sponds to the fact that the magnetic flux associated to 
a magnetic field that is aligned perpendicularly to the 
(xi, a;2)-plane can penetrate this plane with two different 
orientations. 

At zero magnetic field, i3 = 0, the model depends on a 
single parameter [llj . namely, the dimensionful coupling 
constant g with mass dimension 2 — d. In the sense of 
statistical physics, the coupling appears to correspond to 
an irrelevant operator within the perturbative Gaussian 
classification. However, this naive scaling analysis does 
not yield the correct picture for the infrared behavior of 
this model, as we will explain in detail in Sect. I VII From 
the point of view of quantum field theory, the Gross- 
Neveu model provides an ex amp le of a nonperturbatively 
renormalizable field theory [501, and the {i^ipY coupling 
becomes RG relevant near a non-Gaussian fixed point 




III. GROSS-NEVEU SYMMETRIES 

While the construction of local operators that conform 
to Lorentz symmetry can be ensured by taking products 
of operators and performing suitable contractions over 
tensor indices, the transformation properties under dis- 
crete space-time symmetries require more work. 

In Ref. [l^l explicit representations for discrete space- 
time transformations as realized on 4-component spinor 
fields over three dimensional Euclidean space were given 
as follows: 

C : Q = ^ [(1 + 07273 + i(l - O7275] , (11) 
V: Pc = \[{l + C)7i73 + i(l - 07175] , (12) 
r : = i [(1 + vhi + i(l - ry)727o] , (13) 



where C, P, T denote charge conjugation, parity inver- 
sion and time reversal. To each transformation, there 
exists an associated unitary 4x4 matrix, which we de- 
note by Cj, P(; and T^, respectively. Other possible con- 
ventions were given for example in [sij . However, we 
will stick to the definition as displayed in Eqs. ([TT]) . (fT2|) 
and (jlSp . Concerning parity, it is worthwhile to men- 
tion that in three space-time dimensions, parity inver- 
sion is properly defined by (xq, xi, 2:2) f-J- (xq, — xi, X2), 
i.e. only one of the spatial components is reversed. This 
is due to the fact, that in our case, only a single gen- 
erator for rotations exists. The above definition ensures 
that parity inversion is not an element of the connected 
component of rotations containing the identity. As can 
be seen in Eqs. ([Tl]) - (fT3|) . there exists an entire family 
of realizations of discrete transformations depending on 
pure phase variables C ^^'^ V with unit modulus. We 
will simply set ^ = C = 77 = 1 in the following and omit 
the subscripts on the transformation matrices. The re- 
alization of space-time transformations on spinor fields 
reads for charge conjugation 

CV'(x)C-i = ii'Cf, C^{x)C-^ = -{Cmx)f, (14) 

parity 

rip{x)r-^ = PV^i), r^pis:)^^ = 1/^(^)P^ (i5) 

and time-reversal 

T^{x)T-^ = Ti>{x), T^p{x)T-^ = ^{x)T^. (16) 

Here we follow the notation as given in [5l|, and define 

i = (xo, -a:i,X2)^, x = (-xq, xi, ^2)^. (17) 

Note that by virtue of TiT^^ = — i, time-reversal is an 
anti-unitary transformation. We will consider a theory 
to be symmetric under C-, V- and T-transformations, if 
its Lagrangian density obeys 

CC{x)C-^ = C{x), (18) 
VC{x)V-^ = C{x), (19) 
TC{x)T-^ = C{x), (20) 

such that a simultaneous transformation acting on fields 
and coordinates leaves C invariant as a function of space- 
time coordinates. 

The Gross-Ncveu model as defined in Eq. is sym- 
metric under the discrete chiral transformation 

1/' ^ l5i>, "0 '-^ -'075- (21) 

Since the composite operator Tpip transforms into —ipip 
under this discrete chiral transformation, a finite expec- 
tation value {iptjj) 7^ in a given quantum state signals 
the breakdown of chiral symmetry. In the following, we 
denote the discrete symmetry group by = {1,75}. 
However, there exists also a continuous Abelian chiral 
symmetry, generated by 735: 

^^eivj35^^ (22) 
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It is easy to see that the element e'^T^s g U^^(l) com- 
bined with the non-trivial transformation 

i' TaV", ^ -"073 (23) 

leads us back to the transformation Eq. (PT|) . In this 
respect, the discrete does not yield a new symme- 
try of the theory and we will henceforth choose Z2 and 
U'^^(l) to define the chiral symmetries of the three di- 
mensional Gross-Neveu model. It is perhaps worth men- 
tioning that for the free theory, i.e. g = 0, the symmetry 
transformations act independently on all flavor species j, 
J — 1, . . . ,iVf. For finite couplings, however, all flavors 
are subject to a simultaneous chiral transformation, as 
we expect from the symmetry-breaking pattern induced 
by the interaction term. 

To conclude this section, we summarize the properties 
of electromagnetic quantities under discrete space-time 
transformations. The gauge potential A transforms as 

CAf,{x)C-^ = -~A,,ix), (24) 

VA^{x)r-' = A^^ii), (25) 

TAf,ix)T-' = -A^{x), (26) 

where A = (A, -^1,^2)^ and A = {-Ao,Ai,A2 f. 
The field-strength tensor J^^jjy = d^A,^—d,yAf_i accordingly 
obeys 

CJ-^,(x)C-i = -Tf,,{x), (27) 
VT^,,{x)V-' = (28) 
TT^,,ix)T-' = -P^,{x). (29) 

Here, by J-^^ and J-f^^ we denote matrices which result 
from plugging Eqs. (f25|) and ([26|) into the definition of the 
field-strength tensor. The so-called dual 'field-strength' 
Ff^ = \e^u<y^iycr bccomes a pseudo- vector quantity in 
three dimensions. Thus, it behaves as 

CF^{x)C-^ = -F^{x), (30) 

VF^ix)V-' = -F^(i), (31) 

TF^{x)T-' = F^{x). (32) 

The simplest Lorentz, gauge and C-, 7~-invariants 
that can be built from the field strength and its dual are 

and F^ (33) 

Since = ^-^'fi.i' there is only one linearly independent 
invariant. 



IV. CLASSIFICATION OF COMPATIBLE 
OPERATORS 

Having collected the prerequisites to study all oper- 
ators that are in principle compatible symmetry-wise 
with the Lagrangian density Eq. we will give an 
exhaustive classification on the level of fermionic bilin- 
ears in Sect. ITvXl and then proceed to quartic fermionic 
terms including purely magnetically induced operators in 
Sect. HVBl 



A. Bilinear Operators 

Fermionic bilinears are the building blocks of an ac- 
tion for fermionic degrees of freedom. Terms enter- 
ing the quadratic part of the action define the inverse 
bare propagator of the theory. In our case, we need 
to contract both spinor and Lorentz indices to form 
an appropriate scalar quantity. But, due to the pres- 
ence of the gauge background, we can as well use the 
gauge-invariant field-strength tensor and its dual to build 
Lorentz-invariant fermion bilinears. We consider only 
those bilinears which are Lorentz symmetric as a build- 
ing block of the effective action. By use of Eqs. 
([TC]) and Eqs. (PT|) and ([^ we obtain the behavior of a 
given bilinear under discrete space-time and chiral trans- 
formations. The results are collected in Tables HI HIl 
and lllll For turning operators with vector or tensor struc- 
ture into Lorentz-invariant scalars, we contract as indi- 
cated above with F^ or J^^^. By virtue of the canon- 
ical mass dimension of the fermion fields [-0] — 
the mass dimension of a (non-derivative) fermion bi- 
linear is [ipToip] = d-1 with Fs e {1,73,75,735} 
(scalar/pseudo-scalar), Fv € {7^, 737m ' 75 7/^, 7357m (vec- 
tor/axial vector), or Ft G Wf^i,, ^30- fj.„,^5crf^i,, ^35(7 fiu} 
(tensor/pseudotcnsor). Here, 

o'm^' = ^[7m'7i^], /^,'^ = 0,1,2, (34) 

is the set of generators for (Euclidean) Lorentz trans- 
formations. For the contracted terms we obtain 
[^^^(^r^Vi)] = d+1 and [gJ"^y(0Fpi,i/')] = d+1, since 
[Ffj] = [J^i^iy] = I and [q] = Taking into account 

that the space-time integral in a local action contributes 
[J^] = — c? to the total mass dimension of a given local 
operator, these operators naively correspond to irrele- 
vant directions in theory space. The charge q needs to 
be included in this consideration, since we are interested 
in magnetically induced phenomena. We could also per- 
form two contractions of qFfj^ with an appropriate ten- 
sor structure. This would inevitably increase the mass 
dimension by 2 and render this term even more power- 
counting irrelevant compared to mass terms (with mass 
dimension —1) and kinetic operators with mass dimen- 
sion 0. However, from Tables Hl lIIIl we see that none 
of these bilinear operators containing the magnetic field 
are compatible with the symmetries of the Gross-Neveu 
model. They cannot be generated during the RG flow 
in a continuous fashion. Combinations like ■Ff_ii,{il'i7f^^ip) 
and Ff^iTjjjsQj^iJj) conform to space-time symmetries, but 
violate the chiral symmetries of the Gross-Neveu action. 
Chirally invariant combinations containing two fields 
arc not even symmetric with respect to C, 7^ and T- 

To conclude, wc would like to comment that each 
contraction with an appropriate vector bilinear can be 
rewritten as a contraction of J^^j, with a tensor bilinear 
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by means of Eqs. dZ])-®: 

^)'(V'7357/x'0) = 



-^-^M!^ ("0735^^1. -0), (35) 

-^J^^^,^{■ipl5cr^L•^'^P), (36) 

+ ^J^,,,^i-tpricr^,^ilj), (37) 

-^J'tJ.uiip'J^^ip). (38) 



not make use of this full theory space, but only repro- 
duces the Gross-Neveu coupling in the derivative expan- 
sion. This might be accidental or could point to a further 
hidden symmetry. 

We now consider quartic fermionic terms that can be 
built by contracting once with or f/i^^j. Due to the 

identities Eqs. (P5l) - ([55|) . it is sufficient to consider con- 
tractions with qFf^ only. From Tables Hl lIIII we conclude, 
that the following terms are allowed 

ii'i^) (^7357m^m^) = - ^ (^V') (V^(T^.^^. V-) , (42) 



B. Quartic Operators 

Quartic fermionic terms describe the two-body inter- 
action processes in our theory. We first ask for terms 



(43) 



Nf 



(39) 



which obey the symmetries of the original Gross- 
Neveu action and have canonical mass dimension 
[{'iljiTx'>Pi)i'4'j^Ytpj)] = 2(i — 2. In a next step, we an- 
alyze contributions with mass dimension 2d, i.e., quartic 
terms that include a contraction with a mass dimension 
2 object. In principle, one could also move to operators 
with higher mass dimension. But with increasing mass 
dimension our naive expectation is that these operators 
become increasingly irrelevant for the IR behavior of the 
theory, even near a non-Gaussian fixed point. 

Note that terms as captured by Eq. ([59)1 arc composed 
of flavor singlets. Invariant sums of bilinears with off- 
diagonal flavor structure can be brought into this singlet- 
singlet form by an appropriate Fierz transformation. The 
result is in general a sum over several singlet-singlet con- 
tributions with different Lorentz structure. 

Since a discrete chiral symmetry is less restrictive than 
continuous chiral symmetry, the number of allowed quar- 
tic terms appears to be quite large. However, the re- 
quirement for invariance under the continuous U'^^(l) 
symmetry remedies the situation somewhat. The al- 
lowed terms in the absence of a magnetic field are ex- 
hausted by iij^)^, (?^7^V)^i (■00-^i.V')^ and (V'735V')^ 
("07357^0)^? (V'735<''mi'V')^ as well as the invariant com- 
binations (■073'!/')^ + (■075-0)^ > (0737m "0)^ + (V'757m'0)^ 
where we have omitted all flavor indices. Actually 
(V'7357m'0)^ and (■0735O'mi''0)^ can be written as 



(V'7350-Mi'V')^ = ^(V'7M'0)^ 

(V'7357m'/')^ = \{'4>cr,j.uipf 



(40) 
(41) 



since they are simply linear combinations of the Dirac 
algebra basis elements. The remaining six quartic terms 
are mutually independent. They parametrize the 'Gross- 
Neveu theory space' without an external magnetic field 
in the pointlikc limit. It is interesting to note that an 
RG flow calculation in the pointlike limit does actually 



and 



('/'73V') (0757m ^p^) + (V'750)(0'737AffM0)- (44) 

Of course, a U^'^(l) invariant operator with fully con- 
tracted Lorentz indices squared is also compatible with 
the Gross-Neveu symmetries, but has mass dimension 
larger or equal to 2d + 2. One could also consider func- 
tions of the invariant [qF^Y times a bilinear squared. 
This, however, corresponds to operators with arbitrarily 
high mass dimension. In Appendix|5]we show that start- 
ing from the naive Gross-Neveu action Eq. ([T|) , a term of 
the type Eq. (|42p is indeed generated in an infinitesimal 
RG step. In phases with broken chiral symmetry, such 
an operator may give rise to an anomalous magnetic mo- 
ment, as observed in 1271. 



V. FUNCTIONAL RG EQUATION 

In this section we briefly summarize the essentials of 
the (functional) rcnormalization group. In the space of 
couplings gi that parametrize a given action, the RG pro- 
vides a vector fleld /3, summarizing the RG /? functions 
for these couplings (/3)i = /3gi(5i,52, . ■ . ) = dtgi- As the 
full content of a quantum theory can be specified in terms 
of generating functionals for correlation functions or ver- 
tices, we can more generally study the RG behavior of 
a generating functional. Introducing an IR-regulated ef- 
fective average action by inserting a scale-dependent 
regulator into the Gaussian measure of the functional in- 
tegral representing the generating functional of vacuum 
correlators, Zk[J], the RG fiow of this action is deter- 
mined by the Wetterich equation [52| 



{dtRk)\, dt = k 



(45) 

More precisely, Tk [$] is the modified Legendre transform 
rfe[$]-(J,$)-iyfe[J]-i($,i?fc$) (46) 

of the scale-dependent generating functional 

Wk[J]^ log Zk[J] (47) 
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V 


T 








('ipip)(x) 


{Tpip)(x) 


{iiip){x) 


-{i>ip){x) 


(7/)7/))(a;) 


(?/)73V')(2:) 


-{4>'y3Tp){x) 


-{ipyiTp){x) 


-{lpl3i'){x) 


{ip'y3i>){x) 


(7/)73e2''^^35^)(2,) 


(V)75V')(a;) 




{i'l5lp){x) 


-{lpl5lp){x) 


-(V'75V')(a;) 


(7/)75e2''^'^35^)(2,) 


(^735^/')(a;) 




-(i/'735i/')(i) 


-(l/'735l/')(4) 


(V'735i/')(a;) 


(^735 V') (a;) 



Table I. Classification of scalar/pseudo-scalar fermion bilinears {■tl)Tsi}){x) according to their behavior under discrete space-time 
and chiral transformations. 
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T 




U^^(l) 


{■ip^^tlj){x) 




('/'7m 


-(V'7m 


(V'7M'/')(a;) 


{ipf f,ip){x) 


(i/'737MV')(a:) 


-(V'737mV')(2;) 


-(V'737mV')(S) 


(V'737mV')(*) 


-(V;737MV')(a;) 


(^737Aie^'*'^^^i/')(a^) 


(V'757M^)(a;) 


(i/'757M'0)(a;) 


(i/'757MV')(i) 


(i/'757mV')(^) 


(V'757M^)(a;) 


(^'757^6^"^^^^^) (2;) 


(V'7357M'0)(a^) 


-(i/'7357MV')(a;) 


-(i/'7357m 


(V'7357mV')(*) 


-(V'7357M'/')(a;) 


(i/'7357MV')(a:) 



Table II. Classification of vector/axial- vector fermion bilinears {tpTvip){x) according to their behavior under discrete space-time 
and chiral transformations. Here, we have defined 7 — (70,-71,72)"^ and 7 = (—70,71,72)"^- The bilinear "07357/^1/' can be 
shown to be related to -tpa^utp, cf. Eq. ((7|. 
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V 
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U^^(l) 


{il'a^^Tp){x) 


-{■ipa^^Tp){x) 


{Tpaf,„ip){x) 


{ipaf,„ip){x) 


-{ipa^^ip){x) 


['4>CJl,^^l)){x) 


(V'730-M-'V')(a^) 


{i>j3cr^^jp){x) 


-(^735-M-- 


{ipj3^,,„Tp){x) 


{ipj3<yt,^ip)ix) 


(^730-^1^6^"^^=*= ■(/)) (a; 


(V'75crp^V')(2:) 


-{i>To'^^^ip)(x) 


(^755-M'-'/')(i) 


(V'750-M'^V')(a^) 


~{i>T^<J^u.'ip){x) 


(^75cr^i^e^"^^^^?/>)(2; 


ii''y35cr^^Tp){x) 


-{'4>-/35<Jf,^-4)){x) 


-(i/'7355-M,.i/')(i) 


{ip'y35a^^ip){x) 


(V'735cr^^i/')(a;) 


(V'735cr;.,.7/')(a;) 



Table III. Classification of tensor/pseudo-tensor fermion bilinears {'iI}Tt^){x) according to their behavior under discrete space- 
time and chiral transformations. The matrices cr^j^ and a^j^ carry the same sign structure as T^iv and F^j^v in Eq. (|28|) 
and (|29p . respectively. The bilinear ^jj'ysia/j.^ijj can be shown to be equivalent to V'7^'i/'i cf. Eq. ((9|. It also holds that 

^73/50-Mf ~ ^75/37mV'- 



for connected correlators and Tf, is the second func- 
tional derivative with respect to the field representing 
an appropriate collection of field variables for all bosonic 
or fermionic degrees of freedom. The brackets ( • , • ) are 
shorthand for suitable contractions of generalized indices. 
The function Rk denotes a momentum-dependent regu- 
lator that suppresses IR modes below a momentum scale 
k. STr denotes the so-called super-trace operation, which 
simply takes into account the minus sign from the loop 
contributions from the fermionic sector. The solution 
to the Wetterich equation provides for an RG trajec- 
tory in the space of all action functionals, also known 
as theory space. The trajectory interpolates between the 
bare action S\ to be quantized Tk^\ — > Sa and the full 
quantum effective action F = Tk^Q , being the generating 
functional of IPI correlation functions; for reviews with a 
focus on fermionic systems, see |53l456| and particularly 
\Et\ . Parametrizing the effective average action by a 
set of generalized dimensionless couplings gi, the Wet- 
terich equation provides us with the corresponding RG 
flow dtgi = j3g\gi, g2, ■ ■ ■)■ A fixed point gi * is defined 
by 



ft(5i,*,52,*,...) = 0Vz. (48) 



In general, this map of a quantum field theory to a fiow 
in coupling space allows to extract universal physical in- 
formation from the fixed points and from the associated 
manifolds controlling the flow toward the infrared. 

VI. MAGNETIC BETA FUNCTION 

In the following, we will elaborate on how the phe- 
nomenon of magnetic catalysis manifests itself in the IR 
behavior of the beta function for the coupling g. The 
beta function for vanishing gauge background was de- 
rived previously within the formalism of the functional 
rcnormalization group [Tl| . Explicitly, it is given by 

/3, ^ d^g ^ (d - 2)5 - 4 (^^^^) V, /f (0) g' (49) 

with the dimensionless coupling 

9 = k^-^-g- (50) 

In Eq. (|49l). d-y = trlL4 = 4 for the present reducible 
fermion representation, and vj^ = 2''"'"^7r''/^r(d/2). For 
the definition of the threshold function (w) parametriz- 
ing the RG-scheme dependence, see Appendix [Cj Aside 
from the Gaussian fixed point of a nonintcracting theory. 
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Figure 1. The upper part depicts a diagrammatic represen- 
tation of the beta function for the dimensionless couphng g 
(red dot). The regularized fermionic loop contains the full 
B dependence as given by the fermionic propagator coupling 
to the external gauge field A to all orders as shown in the 
lower part. The coupling to A is indicated by a cross. The dt 
derivative acts on the regulator-dependent part of the propa- 
gators and generates diagrams with regulator insertions (for 
details, see Appendix [B|). 



this beta function has also a non-Gaussian, i.e., strong- 
coupling fixed point. At this non-Gaussian fixed point 
the scaling dimension is altered from 0Gauss(£) < at 
the Gaussian one to 9nonGauss(ff) > 0, cf. [ll|. So in 
contrast to a naive scaling analysis, the coupling actu- 
ally corresponds to a relevant direction in theory space at 
this non-trivial zero of the beta function. This fixed-point 
structure is related to the known nonperturbative renor- 
malizability of the Gross-Neveu model. The renormal- 
ization program was performed within a 1 /TVf-expansion 
[soj . and it was argued that the theory is renormaliz- 
able to all orders in this expansion. With the help of 
the renormalization group, the theory was shown to be 
asymptotically safe, i.e., the RG flow contains trajecto- 
ries that eventually hit the non-trivial fixed point with 
finite dimensionless couplings if the scale k (and implic- 



itly the UV cutoff A) is sent to infinity In this 

model, the non-trivial fixed point also acts as a sepa- 
ratrix for different regimes of the fiow. These two dif- 
ferent regimes translate into different phases realized in 
the respective ground states of the quantum field the- 
ory [s^, nil . For values g < g*, the flow is toward a non- 
interacting theory, while for g > g* the coupling diverges 
at some finite IR scale fee, indicating that the discrete 
chiral symmetry is spontaneously broken. Thus, can 
be interpreted as the critical value of g corresponding to 
a chiral quantum phase transition. Whereas g^ is gen- 
erally non-universal, the universal (scheme-independent) 
critical exponents of this quantum phase transition can 
be determined quantitatively with v ery good agreement 
among the various methods [ll|, l59l - [6l| . The order pa- 
rameter for chiral symmetry breaking is the so-called chi- 
ral condensate, i.e., the expectation value of the operator 
Sil^i It can be convenient to introduce an order- 

parameter field a into the functional integral by means 
of a Hubbard-Stratonovich transformation. However, for 
giving a simple renormalization group picture of mag- 
netic catalysis and also for avoiding technicalities, we 
stay within a pointlikc, purely fermionic description. For 



finite flavor number, we expect our arguments to be reli- 
able until the onset of chiral symmetry breaking, where a 
proper description of the built up of a finite expectation 
value of the order parameter and the associated collective 
fluctuations need to be taken into account. 

To capture the effects of a magnetic field on the RG 
flow, we make the following ansatz for the effective aver- 
age action: 



rk[i}j,^j,A]= (51) 



where the coupling g is now to be understood as a func- 
tion of the RG scale k. Inserting this ansatz into the ffow 
equation, we can project onto the flow of the dimension- 
less coupling g by taking appropriate functional deriva- 
tives with respect to the spinor fields. The presence of 
the magnetic field requires a slightly different approach 
in the derivation of the beta function, which is usually 
performed for the corresponding momentum-space quan- 
tities and fields. One problem is that the vector po- 
tential A{x) entering the quadratic part of the action 
renders the fermionic propagator a translationally non- 
invariant function of the space-time coordinates. Tak- 
ing gauge holonomy factors appropriately into account, 
it still becomes possible to study the ansatz in Eq. (|5ip . 
which essentially constitutes the zeroth order in a deriva- 
tive expansion. This implies that the derivative structure 
contained in the effective average action is the same as 
in the classical action. In Appendix [B] we demonstrate 
that Eq. ([?T|) leads to a consistent fiow equation within 
such a derivative expansion. There we also show that the 
scale dependence of a wave function renormalization en- 
tering the action in the form Z^p V'ji^l-^] V'j in the 
present truncation is trivial, i.e., = 1 like in standard 
flows for pointlike fermionic theories with chiral symme- 
try [Hill. 

The beta function /3g = dtg is derived in Appendix IB] 
It incorporates the coupling of the fermions to the exter- 
nal field to all orders in qB, see Fig. [TJ and reads using 
a Callan-Symanzik regulator 



I 2' 26 



26 



(52)" 



Here, we have defined the dimensionless magnetic field 
parameter 



b = 



qB 

fc2 ■ 



(53) 



In the limit 6 — > 0, the correct results for vanishing field 
are recovered, cf. Eq. ([i^ and Appendix The mag- 
netic field defines a new scale, the so-called magnetic 
length I — 1/^/qB. The inverse magnetic length thus 
describes the typical energy and momentum scales as- 
sociated with fermions moving in a magnetic field. For 
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Figure 2. Plot of the beta function Pg (cf. Eq. ((52)) ') of the 
coupling g in the presence of an external magnetic field B in 
arbitrary units. As the RG scale k tends to zero for fixed 
B, the quadratic part of the beta function dominates over 
the linear part due to dimensional flow. The strong-coupling 
fixed point (black dots) is pushed toward the Gaussian fixed- 
point as the scale moves from the high UV to the IR regime 
(indicated by dashed red arrow). This leads to a divergence in 
g at a finite scale associated with chiral symmetry breaking, 
for arbitrary values of the initial coupling at the UV scale A. 
From right to left h takes the values 0.1, 1, 2, 3, 5, 10, 100 
and 1000. 



so-called "quantum anomaly" [16[ for the Lagrangian 

C^Y.'^imA]^,. (54) 

It can be shown that the U(2A'f) symmetry of this 
noninteracting system minimally coupled to a perpen- 
dicular magnetic field is spontaneously broken down to 
U(iVf) X U(iVf) for B 7^ 0. One finds EUl-iil 

^(v:;,^,) ^ -iVf|-. (55) 

However, for vanishing interaction (g = 0), this expecta- 
tion value can not couple back into the bilinear part of 
the theory. So, although the chiral symmetry is broken 
and a finite chiral condensate is dynamically generated, 
it does not act as a mass term for the fermions and the 
spectrum of the theory remains unchanged. 

VII. RG ANALYSIS OF MAGNETICALLY 
INDUCED SPECTRAL GAP 



strong magnetic fields, the fermionic fluctuations that ex- 
perience Landau level quantization are primarily driven 
by the lowest Landau level [n = 0, t = — 1 (cf. Ap- 
pendix [^l]- A Landau level can roughly be attributed an 
extent of I in all spatial directions orthogonal to the mag- 
netic field. The present quantum field theory in d = 3 is 
thus dimensionally reduced via the external field to es- 
sentially a quantum mechanical problem (d — > d—2 = 1). 

In the regime <^ k, the high-momentum quantum 
fluctuations on the scale k arc not affected by the long- 
wavelength physics associated with the low-momentum 
scale From the initial UV scale k = A down to the 
magnetic scale k = we expect the beta function to 
be very similar to its zero field counterpart. However, as 
we evolve the system toward the IR by sending the RG 
scale k to zero, the IR fluctuations in the regime ^ k 
are dominated by magnetic-fleld effects. 

This behavior is depicted in Fig. [2] At high scales 
(6 3> 1), the beta function shows the typical competition 
between a positive linear and negative quadratic contri- 
bution [cf. Eq. (|49| ]. Its shape is unaffected by the pres- 
ence of an external magnetic field. Upon lowering fc, the 
non-Gaussian fixed point starts moving toward the Gaus- 
sian fixed point. In the deep IR limit, the non-Gaussian 
fixed point asymptotically merges with the Gaussian one. 
However, only the negative branch of the beta function 
survives in this limit, rendering g a relevant perturba- 
tion, completely oblivious of the initial UV value of the 
coupling. In this way, the divergence of g at a finite scale 
is catalyzed by the external field for arbitrary values of 
9- 

A complementary approach to this phenomenon is the 



The phenomenon of spontaneous chiral symmetry 
breaking is accompanied by the dynamical generation of 
a mass rhf for the fermionic degrees of freedom which 
drive the symmetry breaking. In a partially bosonizcd 
formulation, which we already alluded to earlier, the con- 
densation of an order-parameter field a introduces a mass 
term into the free part of the Dirac Lagrangian. The on- 
set of symmetry breaking driven by quantum fluctuations 
typically occurs at some critical scale k^. We expect that 
fhf ~ kc- The collective fluctuations of the fermionic 
system encoded in fluctuations of the a field about its 
expectation value will typically diminish this value. In 
cases where a continuous symmetry is broken, the ac- 
companying Nambu-Goldstone modes have the tendency 
to restore the order, and the expectation value of the 
order-parameter field can be pushed to smaller values. 

The beta function Eq. can now be used to analyze 
the field dependence of the dynamically induced fcrmion 
mass. For this, we use the identification rhf = kc, ignor- 
ing possible 0(1) factors arising from fluctuations in the 
symmetry-broken regime. This approximation is equiva- 
lent to a gap-equation approach [62[. The critical scale 
kc can be identified with the scale where the flow of g di- 
verges. Therefore, it is convenient to reformulate the beta 
function for g as a beta function for the inverse coupling 
l/g. From /3i/c,, we can move to the actual trajectory 
\/g{k) by solving the fiow equation with a given initial 
condition g{h) = go and record whenever \/g[kc) = 0. 

The result of this analysis is shown in Fig. [31 The 
initial or bare value go was chosen from three different 
regimes. The flow equation was evaluated for flavor num- 
bers ranging from A^f = 1 to A^f = 100, which at this level 
of the truncation is indiscernible from the large- A^f limit. 
For clarity, only the curves for A'f = 1 and iVf = 100 are 
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Figure 3. B dependence of the dynamically generated fermion 
mass rh{ as determined from the critical scale fee- The lowest 
boundary of each band corresponds to Aff = 1 and the upper 
boundary to A*'f — 100. The curves for intermediate A'f can 
be computed as easily, but are not shown for clarity. Band 
(A) is the weak coupling go <^ g*, band (B) the intermediate 
but sub-critical coupling go 5* , and (C) the strong coupling 
regime, go > 5.. 

shown. The curves for 1 < A'f < 100 fall into the shaded 
regions and do not cross in the B-interval we scanned. 
The precise values for go and A arc listed in Table IIVI In 
the weak coupling regime go <^ 5* we find mf cx B. For 
intermediate, but sub-critical values of the initial cou- 
pling go ^ g* deviations from the linear behavior show 
up. In a sense, in this regime the generation of a fermion 
mass is due to the interplay of the 'quantum anomaly' 
of the free system for i? 7^ mentioned in Sect. IVII and 
the four-point vertex g (in the pointlike approximation). 
While the anomaly generates a finite chiral condensate, 
the four-point vertex allows for coupling this condensate 
back into the single particle spectrum, such that a mass 
gap opens up. In the strong coupling regime g ^ g* the 
fermion mass is generated already for i? = in the usual 
chiral quantum phase transition from a massless phase 
to a phase with massive Dirac fermions controlled by the 
beta function for vanishing field. 



Regime 






(DS) 
% 




(A) go < g. 


0.20 


100 a.u. 


0.206 


100 a.u. 


(B) go < g* 


2.00 


100 a.u. 


2.76 


100 a.u. 


(C) go > g* 


6.50 


100 a.u. 


65.00 


100 a.u. 



Table IV. Initial or bare value of the interaction go and cutoff 
scale A in the large-A*'f RG and Dyson- Schwinger gap equa- 
tion. The scheme dependent critical values are gi^*^"* = vr and 
gi°"' 5.39. 



To assess the validity of the pointlike approximation of 
the four point vertex and the naive identification of the 
dynamically generated fermion mass rhi with the critical 
scale fee we use in this work, we compared the large-A'f 
limit of our flow equation with the large- A'f gap equation. 



Figure 4. Quantitative comparison of the B dependence of 
the dynamically generated fermion mass from the large-A^f 
flow equation (solid blue) and the large-A^f gap equation in 
the weak-coupling regime. For cutoff A and bare value of go, 
see Table HVl 



This equation was derived in e.g. [T^ and can be written 
as 

1 = / —e-'^'^qBt) coth {qBt) . (56) 
27r2 J1/A2 t2 

The integral is regularized in the proper-time formula- 
tion and can be performed by analytic continuation, such 
that a closed but cutoff-dependent expression is obtained 
(see ). The results are shown in Figs. HUH) The val- 
ues of the cutoff scale A and the bare interaction go used 
in the two different schemes (RG and Dyson-Schwinger 
gap equation) are given in Table ITVl In order to compare 
these two approaches, the values of A and go were cho- 
sen, such that the dynamically generated fermion mass 
fhf obtained from the RG flow and the self-consistent 
solution to Eq. ([55)) coincide for qB/A^ = 0.1. Since 
the Gross-Neveu model (or its corresponding universality 
class) only has a single parameter (which is true beyond 
the pointlike approximation), this fixing prescription is 
sufficient. 

We find excellent agreement over a wide range of qB. 
Deviations occur due to cutoff effects when qB / A? ~ 
0.35, as expected for scheme dependencies. Unfortu- 
nately, this quantitative check cannot be carried over to 
finite Af in a fully controlled manner. Especially the 
identification rfif = fee cannot be expected to hold any 
longer, since collective excitations - which arc completely 
damped in the large- A'f limit - will modify the value of 
TOf beyond the critical scale fee that for finite flavor num- 
bers signals only the onset of chiral symmetry breaking. 
We also expect the strength of the transition to be weak- 
ened upon inclusion of order-parameter fluctuations, as 
has been recently observed in 4d quark-meson models 
[6|. 
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Figure 5. Quantitative comparison of the B dependence 
of the dynamicaUy generated fermion mass from the large- 
A*'f flow equation (sohd blue) and the large- A'f gap equation 
(dashed-dotted, magenta) in the intermediate to strong cou- 
pling regime. The coupling is still sub-critical. For cutoff A 
and bare value of go , see Table IIVI 
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Figure 6. Quantitative comparison of the B dependence of the 
dynamically generated fermion mass from the large-A'f flow 
equation (solid, blue) and the large- A^f gap equation (dashed- 
dotted, magenta) in the strong coupling regime. For cutofl' A 
and bare value of go, see Table ITVl 



VIII. CONCLUSIONS AND OUTLOOK 

In this paper, we have studied the phenomenon of mag- 
netic catalysis within the three-dimensional Gross-Neveu 
model. Minimally coupling the Dirac fermions to an ex- 
ternal magnetic field catalyzes chiral symmetry break- 
ing for arbitrary bare values of the fermionic interaction. 
Even in the free system, chiral symmetry breaking occurs 
and can be attributed to a form of 'quantum anomaly'. 
However, only in the interacting case does the symme- 
try breaking lead to a mass generation for the fermionic 
degrees of freedom. Using the functional RG flow equa- 
tion for the generating functional of IPI vertices, we have 
obtained the beta function for the coupling parametriz- 
ing the fermionic interaction in a pointlikc approxima- 



tion. An analysis of the IR behavior of the beta function 
yielded a clear picture of magnetic catalysis in the lan- 
guage of the renormalization group. The strong-coupling 
fixed point is pushed toward the Gaussian fixed point. 
But in the deep IR limit, we find that only the negative 
branch of the beta function remains, and subsequently 
the coupling always fiows into the strong coupling regime 
and eventually signals the breakdown of chiral symmetry 
in a divergence at a finite critical scale. 

This fixed-point picture of the approach to criticality 
in fact occurs in many fermionic systems [57| . For in- 
stance, also a gauge-field background giving rise to a non- 
vanishing Polyakov loop pushes the strong-coupling fixed 
point toward the Gaussian fixed point j67| . The unique 
feature of a magnetic background is the occurrence and 
subsequent IR dominance of the lowest Landau level act- 
ing as a symmetry-breaking catalyzer [l^ [13] . 

By a symmetry analysis we provide insight into the 
structure of the Gross-Neveu theory space and more im- 
portantly how additional magnetically induced operators 
might affect chiral symmetry breaking in external fields. 
We find that a fermionic tensor-scalar contribution is gen- 
erated for finite flavor number, when starting from the 
naive Gross-Neveu action with only scalar-scalar interac- 
tions. It is subject to current investigation to study quan- 
titatively the impact of these induced operators on the 
value of the chiral condensate and the dynamically gen- 
erated fermion mass in chirally symmetric quantum fleld 
theories. The occurrence of these magnetically-induced 
operators is not particular to the Gross-Neveu model but 
certainly a generic feature of any fermionic system. In 
fact, this is yet another example for the proliferation of 
operators in the presence of further Lorcntz tensors as 
is familiar from the classification of operators, e.g., at fi- 
nite temperature [13, [sl] . It should be emphasized that 
mean-field studies typically ignore contribution of such 
terms. 

Within our truncation, there is no critical magnetic 
field, i.e., symmetry breaking occurs for all B. In prin- 
ciple, however, once further interaction channels have 
opened up, there might be the possibility to compen- 
sate the 'driving force' from the scalar-scalar channel, 
such that for certain regimes, a critical field Be might 
be necessary to enter the broken phase. This would also 
imply that the interacting theory is not connected to the 
free theory in an analytic way as our present analysis 
suggests, and new scaling laws could be observed. 

Also, we neglected the back-reaction of collective neu- 
tral excitations onto the fermionic system. For small 
flavor numbers, we expect sizable quantitative modiflca- 
tions from the fluctuations of the order-parameter field. 
For infinite fiavor numbers, we recover the exact results of 
the large- A'^f theory. The derivation of a consistent set of 
renormalization group equations for the case of charged 
fermionic and neutral collective excitations is a work in 
progress and will be presented elsewhere. 

Concerning condensed matter applications, precise 
quantitative control over the field dependence of the 
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fermionic single-particle gap might provide understand- 
ing of the interplay between magnetic catalysis and 
interaction-driven anomalous Hall plateaus in quantum 
Hall measurements. 
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Appendix A: Spectrum of Dirac Operator in an 
External Magnetic Field 

Here we briefly summarize essentials of the spectrum 
of the square of the operator i.^[y^], such that the trace 
appearing in the fermion loop (cf. Fig. [1} can be conve- 
niently expressed in the corresponding basis. With the 
conventions in Sec. |TT1 the squared Dirac operator can be 
written as 



D[Afti 



{D{J^ To ® To + qB To ® T3) , (Al) 



where is the covariant Laplacian. The spectrum 

of the squared Dirac operator is thus found to be 

{vl + 4 + TgS I po e n e IN, T = ±1} . (A2) 

Here, = \qB\(2n -I- 1) corresponds to Landau level 
energies and n is the associated Landau level index. Ac- 
counting for the density of states ^ for each Landau 
level, the trace operation can be decomposed as 



27r 



tr(.), (A3) 



where denotes space-time volume and tr is the trace 
over spinor indices. 



Appendix B: Derivation of the beta Function 

It is convenient to derive the flow equations in position 
space rather than momentum space due to the presence 
of the vector potential A(x). Strictly speaking, the prop- 
agator is in general no longer a translationally invari- 
ant function of space-time coordinates. It can, however, 
be decomposed into a holonomy factor times a transla- 
tionally invariant factor [65j . allowing in principle for a 
momentum-space formulation. Still, a formulation in po- 
sition space is more direct in our purely pointlike approx- 



imation. The flow equation Eq. (|45p can be expanded as 

-1 



V k 



isTr 



[dtRk) 



(Bl) 



= -iTra. InG- - ^Tra, ^ ^-^{0,^ 



1=1 



where the dt derivative only acts on the regulator depen- 
dent part. Here we have decomposed the fluctuation ker- 



nel r 



(2) 



-f-Fl , where T)^' is the fleld-independent 



^(2) 



part and V), denotes the fleld-dependent vertex part. 
We also defined the scale-dependent regularized propa- 
gator 



.(2) 



Rk 



We choose a chirally invariant regulator 



Rk — 



fc2 



(B2) 



(B3) 



with an as yet unspecified shape function r{x). The fluc- 
tuation kernel can conveniently be derived in Nambu rep- 
resentation for the spinor fields (see e.g. [11]). We thus 
obtain 



Gk = 



G+ 
Gl 



with Gl 



[Z^,{iIp[A) + Rk]-\ Gl = [G+]^ and 



f(2) ^ [H^ ~Fl^ 



(B4) 



(B5) 



While we have omitted flavor indices for the propagator 
part, which is diagonal in flavor space, the flavor struc- 
ture of the vertex part is non-trivial. The matrix blocks 



are given by 






9 

Nf 


-F'^- = 


^\ 

Nf 1 




Nf \ 




9 

Nf 



(B6) 
(B7) 
(B8) 
(B9) 



We deflne the position space representation of the oper- 
ator i.^[.4] as 

(BIO) 

Accordingly, we proceed for the regulator part Rk in- 
serted into the Gaussian measure of the generating func- 
tional of Green's functions: 



{x\Rk\x') = Rk{x,x') = Z.^, 



ip[A]r 



{x,x'). 
(Bll) 
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In this notation, the regularized propagator is given by 



G+(x,x') = ^P[A] (^1 + jj (,,^'). 

(B12) 

However, the wave function renormalization is = 1 in 
this pointUke truncation. This can be seen from the tad- 
pole diagram giving rise to self-energy corrections. The 

I 



regularized tadpole diagram is proportional to 

dttT:G+{x,x). (B13) 



Since it can be shown that the above mentioned holon- 
omy factor becomes trivial for x' = x, the remaining con- 
tribution to the trace comes from the translation invari- 
ant part and thus the tadpole carries no net momentum 
at all, like in the B = {) case. 

The vertex part of the fluctuation matrix becomes in 
position space representation 



-A, XT 

Nf 



tpf {x)ipj{x)S{x + x') 

Nf 



= -|- ^'ipi{x)ipj{x)+'^^^^{iljiilji){x)5ijl4js{x~x'), 
Hij{x,x') = -j^ipi{x)ipj{x)6{x + x'). 



(B14) 
(B15) 
(B16) 
(B17) 



The projection onto the beta function of the dimension- 
less coupling g = k'^^^g is facilitated by collecting all con- 
tributions on the right-hand side of Eq. (jBip , which are 
quartic in the fermionic fields. We find two contributions, 
iTr^t {G^FG^F} and ^Trdt {CIHG^H}. These sub- 



sequently will be evaluated for spatially constant spinor 
fields. Using the identity 



dtGt 



and after some algebra, we arrive at 



(B18) 



1 - 2 ^' /■ 

-Trd,{GtFGtF + GtHG-H] = —f ^ 



mA] 1 



fc2 



-3 



dtRk {x,x)ipi{iljjtpj) 



fc2 



-3 



dtRk 1 (a;, a;) X! 

(B19) 



where tr denotes the trace over spinor indices. The first 
term on the right-hand side in the first line of Eq. (|B19[) 
contains not only a contribution to the flow of g, but 
also a term ^ (V'iO'^i/^^t/V'i)(V'j'0j) is generated in an 
infinitesimal RG step. Nevertheless, in the large- A'f limit, 
this contribution drops out from the fermionic loop. It 
is also important to note that the spinors ipi and "04 a-rc 
to be contracted with the term in round brackets. To 
see this more clearly, it is suitable to employ a Laplace 



transformation, 

/"OO 

fix)^ / ds/(s)e--^ (B20) 
Jo 

with X ~ {ip[A])'^ / k'^ , in order to handle the appear- 
ing traces. Since the spectrum of (l^[^])2 is known (cf. 
Eq. (|A2[) ). we perform the trace with help of Eq. (|A3p 
and obtain 



qB 
2tt 



A^/tFs sinh 



\ k^ 



cosh I I To (g) To -I- smh ( -p^- ) To ® T3 



fc2 



k^ 



(B21) 
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where the prime indicates that the trace over spinor in- 
dices is not included. The last term in Eq. (|B2ip makes 
the generation of a magnetically induced operator ex- 
plicit. However, neglecting this contribution, the flow 
equation for g can be derived in closed form by specify- 
ing the shape function r{x) as 



r(x) 



a; + 1 



(B22) 



i.e., by using the Callan-Symanzik regulator. The 
Laplace transform is known in this case and reads 



Appendix C: Threshold Functions 

The regulator dependence of the flow equation is en- 
coded in dimensionless threshold functions which arise 
from IPI diagrams describing fermionic quantum fluctu- 
ations. In this work we used a Callan-Symanzik regu- 
lator. In the case of vanishing magnetic field, the beta 
function for the coupling g is given by 



dtg = (d - 2)g - 4 



^Nfd-y - 2 



vdl[{0)g\ (CI) 



with the threshold function 



fis) 



(B23) 



Upon performing the Laplace integral, we finally arrive 
at the rather simple result 



dtg = {d-2)g- 



1 / Nfd-y - 2 

16^ V 



3 J_ 

2' 26 



-2b 



(B24) 

where b = ^ denotes the dimensionless external field 
and y) is the Hurwitz zeta function. 



/f (w) ^ dt / dyy^ 
Jo 



(l + r(2/))2+c.' 



(C2) 



see e.g. [l^l53| for a derivation. For the Callan-Symanzik 
regulator, the threshold function evaluates to lf{0) = \. 
We checked that as & — >■ the beta function Eq. ([5^ has 
the correct limit and the results for vanishing field are 
recovered. This is easily seen by noting that 



lim 

6^0 



= 4. 



(C3) 
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